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Abstract 

This paper derives high-SNR asymptotic average error rates over fading channels by relating them 
to the outage probability, under mild assumptions. The analysis is based on the Tauberian theorem 
for Laplace-Stieltjes transforms which is grounded on the notion of regular variation, and applies to a 
O | wider range of channel distributions than existing approaches. The theory of regular variation is argued 

to be the proper mathematical framework for finding sufficient and necessary conditions for outage 
events to dominate high-SNR error rate performance. It is proved that the diversity order being d and 
the cumulative distribution function (CDF) of the channel power gain having variation exponent d at 

OS '. 

imply each other, provided that the instantaneous error rate is upper-bounded by an exponential function 
of the instantaneous SNR. High-SNR asymptotic average error rates are derived for specific instantaneous 

o 

error rates. Compared to existing approaches in the literature, the asymptotic expressions are related 
to the channel distribution in a much simpler manner herein, and related with outage more intuitively. 
The high-SNR asymptotic error rate is also characterized under diversity combining schemes with the 
channel power gain of each branch having a regularly varying CDF. Numerical results are shown to 



corroborate our theoretical analysis. 
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I. Introduction 

The concept of diversity plays a critical role in communications over fading channels by 
quantifying the decrease of the average error rate as a function of the average SNR [1]. In [2], 
the authors show analytically that if the channel power gain has a probability density function 
(PDF) which behaves asymptotically like a multiple of its (d— l)-th power near 0, the diversity 
order is d. The approach employed in [2] is also utilized in [3] to analyze the SNR shift in more 
complicated systems with the same diversity order, and in [4] to derive the optimal diversity and 
multiplexing trade-off for generalized fading channels. A similar definition for diversity order is 
also established in [5], and used to determine the diversity order associated with some lattice- 
based MIMO detection schemes [6], [7] with equal numbers of transmit and receive antennas. 

Regular variation is a concept in real analysis which describes functions exhibiting power 
law behavior asymptotically near zero, or infinity, and is applied to several different areas 
including probability theory [8], [9]. The Tauberian theorem for Laplace-Stieltjes transforms 
[9, p. 37] asserts that if a function with the non-negative support is regularly varying at the origin 
(infinity), then then its Laplace-Stieltjes transform must be regularly varying at infinity (origin). 
The applications of this in communications and networking are primarily seen in asymptotic 
queueing analysis [10], [11]. In addition, fading channel distributions with properties related 
to regular variation have been studied with emphasis on scaling properties of ergodic channel 
capacity, under several communication scenarios involving channels with heavy tail behavior in 
[12], where capacity scaling of systems rather than diversity analysis is considered. 

To the best of our knowledge, diversity -related performance analysis over fading channels has 
never been addressed using the Tauberian theorem together with the theory of regular variation. 
In this paper, we perform analysis of diversity, and more generally the high-SNR asymptotic 
error rate for a wide range of channel distributions and modulation types under mild assumptions 
on the instantaneous error rate and channel distribution functions. We prove that the diversity 
order is d if and only if the CDF of the channel power gain has variation exponent d at the 
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origin, provided that the instantaneous error rate is upper bounded by an exponential function 
of the instantaneous SNR. Furthermore, we derive more explicit closed-form expressions of the 
asymptotic error rates for the special cases of practical instantaneous error rates that capture 
systems including, but not limited to, M-PSK and square M-QAM. We also establish a unified 
approach to characterize asymptotic error rate performance for diversity combining schemes, 
with the channel power gain of each diversity branch having a regularly varying CDF. The 
results in this paper establish a mathematical framework for determining the conditions under 
which the outage event dominates the error rate performance. Compared to existing approaches, 
our asymptotic average error rate characterization applies to a more general set of channel 
distributions, is related to the channel CDF in a simple manner, is intuitively linked with outage, 
and draws from the well-established mathematical theory of regular variation. Compared to the 
conference version [13], this paper contains complete proofs, and more general assumptions 
on the instantaneous error rate. Furthermore, expanded numerical results not included in [13] 
illustrate improved accuracy compared to [2]. 

In Section II, the mathematical preliminaries regarding regular variation and Tauberian the- 
orem are presented. In Section III we establish the channel and system model through the 
assumptions on the instantaneous error rate function and the channel distribution. Section IV 
establishes a definition of diversity order based on regular variation, which, under general 
conditions, is proved to be equivalent to the most general definition used in the literature. In 
Section V we go beyond diversity analysis and establish the asymptotic equivalence between 
the average error rate and outage. The asymptotic average error rate expressions for diversity 
combining schemes are derived in Section VI, and Section VII concludes the paper. 

II. Mathematical Preliminaries 

Before the mathematical preliminaries, we have a few remarks about notations. Asymptotic 
equality Hi(x) ~ H 2 (x) as x — > a means that \im x ^ a H 1 (x) / H 2 (x) = 1, and Hi(x) = o(H 2 (x)) 
as x — > a means that \im x ^. a Hi(x) / H 2 (x) = 0. E[-] denotes the expectation over the channel 
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distribution, log(-) denotes natural logarithm, and := denote equivalence by definition. Finally, 
7 (a,6) := j^t^e-'dt, T(a) := 7(0,00), Q(x) = ^n)' 1 / 2 exp{-u 2 /2)du, and #„(•) 
denotes modified Bessel function of the second kind of order v. 

We next provide a brief sketch of the mathematical concepts (slow/regular/rapid variation) and 
theorems (a representation theorem and the Tauberian theorem), as well as related propositions, 
which are drawn from standard mathematical references such as [8], [9]. 

Definition 1: A real valued function H(x): R + — > R + is called a Karamata function (at 00 
or 0), if lim^oo H (tx)/H (x) or lirn^o H(tx)/H(x) exists for t > 0. Specifically, the limit 
must be in the form of t rn , where m G [—00, 00] is called the variation exponent of H(x). 

Throughout the paper we interpret t~°° = and t°° — 00 for t > 1, and vice versa when 
t < 1. Specifically, H(x) is termed slowly/regularly/rapidly varying if m = 0, < \m\ < 00, 
and I to I = 00, respectively. By definition, G(x) = H(x~ 1 ) is slowly /regularly /rapidly varying at 
if and only if H(x) is slowly/regularly/rapidly varying at 00. For example, with < m < 00, 
x m is regularly varying with exponent m at both and 00; x m e~ x is regularly varying with 
exponent m at and rapidly varying with exponent —00 at 00; log(l + x rn ) is regularly varying 
with exponent m at and slowly varying at 00. 

For H(x) being slowly varying at 00, there is a representation theorem [9, Theorem 1.3.1, 
p. 12] with the proof available therein. More generally, this theorem can be extended to handle 
the cases of regular/rapid variation, as described briefly in [9, p. 21] and summarized in [14]. We 
express it in our notations as the following, which will be useful in establishing the relations 
among the different definitions of diversity order in Section IV. 

Theorem 1: H(x) has variation exponent m at 00, i.e. Hindoo H(tx)/H(x) = t m for m e 
[—00, 00], if and only if H(x) = c(x) exp { f*(e(u)/u)du} for some a > 0, c(x) — > c e (0, 00), 
and e(x) — > m, as x — > 00. 

It is well-known that the MGF of the channel random variable is central in average perfor- 
mance over fading channels [15]. Instead of the MGF, the Laplace-Stieltjes transform can also 
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be considered, and the theory behind this transform can be brought to bear. The primary tool 
used in this paper is the Tauberian theorem for Laplace-Stieltjes transforms, which relates the 
asymptotic properties of a function and of its Laplace-Stieltjes transform, with both properties 
characterized by slow/regular/rapid variation. This function will often (but not always) be the 
CDF of the channel power gain in the sequel. We summarize [8, Theorem 1, p.443] and [8, 
Theorem 2, p.445] into the following, with the proofs available therein. 

Theorem 2: If a function H(x) > defined on x > has a Laplace-Stieltjes transform 
£( s ) = Jo°° e ~ sx dH(x) for s > 0, then for \m\ G [0,oo], H(x) having variation exponent m at 
oo (or 0) and C(s) having variation exponent — m at (or oo) imply each other. In addition, for 
\m\ < oo and l(x) being slowly varying at (or oo), the relations H(x) ~ x m l(x) as x — > (or 
x — > oo) and C(s) ~ T(m + l)s _m Z(s _1 ) as s — > oo (or s — > 0) imply each other. 

We also have the following proposition which follows from [9, p. 27], and will be used to 
establish the equivalence of the CDF- and PDF-based assumptions on the channel. 

Proposition 1: If H(x) is differentiable with h(x) = dH(x)/dx, and h(x) = mx m ~ l l(x) for 
some function l(x) slowly varying at 0, then H(x) ~ x m l(x) as x — > 0. 
The next proposition will be useful to link the average error rate and outage. 

Proposition 2: If H(x) is slowly varying at oo (or 0), and bounded between two constants 
hi and h 2 (where < h\ < h 2 < oo) for sufficiently large (small) x, then it converges to a 
constant in [hi, h 2 ] as x — > oo (or x — > 0). 

Proof: See Appendix A. ■ 

III. Channel and System Model 

We consider average performance for systems over fading channels where the channel can 
be captured by an instantaneous SNR random variable. These include SISO systems, SIMO 
systems with diversity combining, and MISO systems with beamforming. To facilitate subsequent 
derivations, we factor the received instantaneous SNR into channel-independent and channel- 
dependent components by expressing it as a product pz, where the average SNR (per symbol) p 
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is deterministic, and z is the channel power gain random variable having CDF F(z), resulting 
from both the channel and the system setup (e.g. beamforming or diversity combining). The 
instantaneous error rate P e (pz), as a function of pz, is determined by the modulation type 
together with the noise distribution and represents bit, or symbol error rate. 
In this paper, we focus on analyzing the average error rate 



for large p. We make one of the following assumptions on the instantaneous error rate in our 
subsequent derivations 

. ASla P c (pz) < j3e~ apz V p, where a, ft e (0, oo) are constants. In other words, P e (pz) is 

upper-bounded by an exponential function of pz; 
• ASlb P c (pz) = (5e~ apz , where a, /3 e (0, oo) are constants; 

. ASlc P e (pz) = #2(0) exp(— pz/gi(9))d9, where 9 t e (0, 7r), and the corresponding 
gi(9) and g2{9) are both finite non-negative. In other words, P c (pz) is a positive mixture 
of decreasing exponential functions of pz. 
It can be easily shown that ASlb =>- ASlc =>- ASla. We will have results of differing generality 
corresponding to these assumptions. ASla holds for all modulation types with the additive noise 
being Gaussian. ASlb is a special case of ASla typically applying to non-coherent modulations, 
e.g. P c (pz) = (l/2)e~ pz is the BER for DPSK. ASlc applies to many practical settings with 
appropriate choices of 9 t , g x {ff) and g 2 {9). Specifically, the SER of M-PSK [15, (8.22)] fits ASlc 
with0 t = (l-M- 1 )^, gi{9) = sin 2 9/ sin^M^Tr) and g 2 {9) = tt" 1 ; by making use of Q a (x) = 
7r_1 Jo 7/2a exp(-x 2 / '(2 sin 2 9))d9 for a = 1,2 [15, (4.2) and (4.9)], it is possible to express the 
BER P c (p2) = g( v / 2pi) for BPSK and the SER P c (pz) = 4(1 - M~ 1 I 2 )Q{ ^3pz/(M - 1)) - 
4(1 - M- 1 / 2 ) 2 g 2 ( v /3pz/(M- 1)) [15, (8.10)] for square M-QAM in the form of ASlc. For 
an arbitrary two-dimensional signal constellation with polygon-shaped decision regions over 
AWGN, P c {pz) can be treated as linear combination of the error rates associated with individual 




(1) 
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constellation points given by [15, (5.71)], which is also in the form of ASlc. 

We make one of the following assumptions on the channel in terms of the CDF F(z), which 
is assumed not to depend on the average SNR p: 

• AS2a \\m z ^Q F (t z) / F(z) = r d , where r, d > 0. In other words, F(z) is a Karamata 
function of z at with variation exponent d e (0, oo]; 

. AS2b Same as AS2a with d G (0, oo), in which case one can write F(z) = z d l(z) with 
l(z) slowly varying at 0; 

. AS2c The PDF f(z) := dF(z)/dz exists and f(z) = dz^H^z), where d E (0, oo) and 
li(z) is slowly varying at 0. 
Based on Proposition 1, AS2b is implied by AS2c. Therefore the assumptions, as listed, get 
stronger, i.e., AS2c =>- AS2b =>- AS2a. The difference between AS2a and AS2b is the allowance 
of d = oo, which is equivalent to the rapid variation of F(z) at the origin. This holds, for example, 
for log-normal shadowing. When d = oo is ruled out, we have AS2b, which will be seen to 
offer sharper results than offered by AS2a. When the PDF of the channel power gain exists and 
is in a simpler form than the CDF, AS2c can be employed. This assumption is similar to (but 
more general than) that of [2], which is also based on the PDF. 

As mentioned, it can be verified that log-normal shadowing for which \ogz is Gaussian has 
a CDF F(z) satisfying AS2a with d = oo. Also, it is easy to see that Rayleigh, Nakagami- 
K (Ricean) and Nakagami-g (Hoyt) fading channels satisfy all three assumptions with d = 
1. Furthermore, Nakagami-m fading with f(z) = (T(m))~ 1 m m z m ~ 1 exp(— mz) and F(z) = 
^(m,mz) can be verified to satisfy all three assumptions with d = m; Weibull fading defined 
by f(z) = kz k ^ exp(—z k ) or by F(z) = 1 — exp(— z k ) has d = k. Consider also generalized- K 
fading which models composite multi-path fading and shadowing [16]. For such fading with 
f(z) = 2(kmY k+m '>/ 2 z^ +m - 2 ^ 2 K k _ m [2(kmz) 1 / 2 }/(T(m)T(k)) [16, eqn. (2)], m, k > 0, based 
on Kq{x) ~ — log(x/2) — 7 cm (7 cm 0.5772 is the Euler-Mascheroni constant) and K v (x) ~ 
r(|z/|)(2x _1 )l I/ l/2 for v 7^ near x — 0, it can be verified that AS2c holds with d = min(m, k). 
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IV. Three Definitions of Diversity Order and Their Relations 



Consider the following definitions of diversity order: 



P c (p) = (Gp)- d + o(p- d ) 



(2a) 




(2b) 



= -d 



(2c) 



p-s-oo log p 



where G,t G (0, oo) are constants. As mentioned in Section I, (2a) and (2c) are adopted in 
existing literature (e.g. in [2] and [5] respectively), while (2b) is our preferred novel definition 
indicating that P c (-) is regularly varying at oo with exponent — d when d < oo, and rapidly 
varying with exponent — oo when d = oo. Note that the definition (2a) requires d G (0, oo), 
whereas (2b) and (2c) allows for d = oo as well. We interpret the rapid variation (d = oo) 
in (2b) or (2c) as an average error rate which decays faster than p~ m , as p — > oo for any 
m G (0, oo). In other words, d = oo can be interpreted as the P e (p) versus p plot on a log scale 
does not become a straight line, but "curves down" as the average SNR increases. 

There has been some previous work on relating the channel PDF near the origin to asymptotic 
error rates. The seminal work in [2] shows that if f(z) ~ az m as z — > 0, then (2a) holds with 
diversity order d = m + 1 and array gain G = a((2 d - 1 aT(d + l/2))/(/9Vjrd))~ 1/d - A relation 
between average error rate and outage probability is also established by quantifying how the 
respective diversity orders and array gains are related. The approach in [2] assumes (i) the 
existence of the PDF and its Maclaurin expansion near the origin, and is not naturally linked 
with outage; (ii) the instantaneous error rate is in the form of a Q function (a special case of 
assumption ASlc); (iii) the average error rate is of the form in (2a) with d < oo. 

As an example, it can be verified that f(z) ~ 2m 2m z m_1 (-log(m^/i) — "y cm )/T 2 (m) and 
lim^o K z ) = 00 f° r generalized- fading with m = k, and f(z) does not admit a Maclaurin 
expansion, thus the approach in [2] is not applicable to this kind of channel. We believe that the 
effective mathematical framework for studying high-SNR asymptotic error rates is the Tauberian 
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theorem, with a distinct advantage of characterizing when outage events dominate error rate 
performance by linking the CDF F(-) with P e (p) at high average SNR. Our approach for high- 
SNR analysis enables less restrictive assumptions on both the instantaneous error rate and the 
channel distribution compared to [2], as we will highlight in the sequel. 

The following result establishes that, under general conditions, diversity order defined through 
(2c), is equivalent to our definition (2b) when the limit in (2b) exists. 

Proposition 3: If P c (p) is a Karamata function, (2b) is equivalent to (2c) for < d < oo. 
Moreover, definition (2a) is a special case of (2b) and (2c) (for < d < oo). 

Proof: See Appendix B. ■ 

Since all average error rate expressions in existing literature are Karamata functions, the novel 
definition in (2b) is equivalent to the general one in (2c), and subsumes (2a). 

V. Asymptotic Analysis based on the Tauberian Theorem 

In this section, we prove the asymptotic equivalence between the high-SNR average error 
rate and the channel distribution at deep fading (i.e. outage with small threshold), through the 
Tauberian theorem. In Section V-A, we first establish a necessary and sufficient condition for 
P e (p) = J °° P c (pz)dF(z) to exhibit diversity order of d for P c (pz) satisfying ASla. Furthermore, 
the relation between the asymptotic expressions of P c (p) and F(z) is also characterized through 
regular variation for finite d. In Section V-B, we derive the high-SNR asymptotic average error 
rates for the instantaneous error rate satisfying ASlb or ASlc, together with F{z) satisfying 
AS2b or AS2c. These results reveal a convenient and effective way to characterize asymptotic 
error rate performance. Although the relation between error rate and outage is addressed in 
several ways like in [2] and [17], we present brand new way to relate average error rate to 
outage probability by expressing the asymptotic average error rate in terms of a scaled version 
of outage probability. 
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A. Exponentially Bounded Instantaneous Error Rate 

For the general case with P e (pz) satisfying ASla, we have the following theorem which links 
the notion of regular/rapid variation with diversity order, and allows for d = oo. 

Theorem 3: For P c (pz) satisfying ASla, P c (p) exhibits a diversity order of d G (0, oo] if the 
CDF F(z) of the channel power gain satisfies AS2a. The converse holds if F(z) is a Karamata 
function (at 0). 

Proof: See Appendix C. ■ 
Theorem 3 fundamentally characterizes, with sufficient and necessary conditions, the diversity 
order in terms of the CDF (outage) for small arguments. This is unlike [2] which only provides 
sufficient conditions on the channel distribution to achieve a certain diversity order. We now 
have the following theorem which, unlike Theorem 3, rules out d = oo and assumes d < oo, 
but in return provides stronger results about how the asymptotic expressions of P e (p) and F(z) 
are related, thereby characterizing when the outage event dominates the error rate performance. 

Theorem 4: For P c (pz) satisfying ASla, if the CDF of the channel power gain satisfies 
either AS2b or AS2c, P e (p) ~ Cip _d Z(p _1 ) = c 1 F{p~ 1 ) as p — > oo where c x G (0, oo) is 
a constant. Conversely, if P e (p) = p~ d r(p) as p — > oo where r(p) is slowly varying at oo, 
F(z) ~ c 2 z d r(z~ 1 ) = c 2 P c (-2~ 1 ) as z — > where c 2 G (0, oo) is a constant, assuming that F(z) 
is a Karamata function (at 0). 

Proof: See Appendix D. ■ 
Theorem 4 goes beyond characterizing diversity order and points out the asymptotic propor- 
tionality of P c (p) with F(p' 1 ) as p — > oo, for general P e (p^) satisfying ASla. This naturally 
establishes sufficient and necessary conditions on the asymptotic equivalence between average 
error rate and outage. An implication of Theorem 4 is that for two different communication 
systems over the same channel, there always exists a constant SNR offset between their error 
rate performance at sufficiently high average SNR, as long as their instantaneous error rates are 
both exponentially bounded. 
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B. Specific Instantaneous Error Rates 

We have already established the regular variation of the CDF of the channel power gain at 
as a necessary and sufficient condition for a specific diversity order, under general modulation 
types with P c (pz) satisfying ASla, using the Tauberian theorem. We now offer the sharper results 
when P e (pz) satisfies ASlb or ASlc, with the channel distribution satisfying either AS2b or 
AS2c. The following results make stronger assumptions about the instantaneous error rate (ASlb, 
ASlc) but offer closed-form expressions for constants c\ and c 2 in Theorem 4. 

For P e (pz) satisfying ASlb, the asymptotic expression of the average error rate follows 
directly from Tauberian theorem (Theorem 2) since the average error rate is in the form of a 
Laplace-Stieltjes transform. Given AS2b or AS2c, the asymptotic average error rate is given by 

POO 

PcAsib(p) := P / e- apz dF(z) ~ f3T(d + l)(ap)- d Z( a - 1 p" 1 ) = f3T(d + l)F(c*- V 1 ) (3) 
Jo 

as p — > oo, where the asymptotic equality holds due to Tauberian theorem (Theorem 2) with 
F(z) corresponding to H(x), and ap corresponding to s. The last equality in (3) holds since 
F(z) = z d l(z). Equation (3) shows the high-SNR asymptotic average error rate in terms of 
the CDF F(-). Conversely, given the average error rate P e ,Asib(p) with diversity order d, the 
asymptotic CDF of the channel power gain with small argument can be obtained through the 
substitution p = oT x z~ x , based on Tauberian theorem (Theorem 2). 

For P c (pz) satisfying ASlc, the asymptotic expression of the average error rate 

P e ,Asic(p) := J" £g 2 (e)exp (~^y) dddF ( z ) = £ J" 92$) exp dF(z)d6 

(4) 

can be obtained by evaluating lim^oo P e ,ASic(p)/ F{p~ X )- Note that in (4) we exchange the order 
of integrations due to the positiveness of the integrand and the finiteness of the integral [18, p. 457, 

C. 9]. Like in (3), due to Tauberian theorem (Theorem 2) we have G(p) = J °° exp(— pz)dF(z) ~ 
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T(d + l)F(p 1 ) as p — > oo and has variation exponent — d at oo, then 

Um P^ASlcip) = ^ G(P) . Um Pe^lcip) = r( 1} ^ 1 {e)Q f P_\ dQ _ 

p^oc F( P ' 1 ) p^ooF^p' 1 ) p^oo G(p) v ; P ^ooG(p)J ww \gi(9) J 

r(d + 1) T um g2{e) G{p ' 9 ( l ^ )) de = r(d + 1) f\ 2 {9) 9 i{9)d9 

Jo p^°° G IP) Jo 

(5) 

where we changed the order of limit and integral in the third equality based on the uniform 
convergence condition following from [9, Theorem 1.5.2], and consequently we obtain 

Pe,Asic(p) ~ r(d + 1) ^jT * g 2 (0)gf(9)d0^ Ftp" 1 ) (6) 

as p — )■ oo. Like for (3), the converse of the result in (6) (from average error rate to asymptotic 
CDF) is obtainable through the substitution p = z~ x . For many practical cases, the closed-form 
expression of the integral J" 4 g2(0)gf(6)d9 as a function of d can be worked out without an 
integral. For example, by knowing the instantaneous error rates discussed in Section III, we can 
obtain P c ,bpsk(p) ~ T(d + l/2)F(p- 1 )/(2 v ^F) as p -> oo for BPSK. 

In summary, for the two cases analyzed above, the asymptotic error rate expressions in (3) and 
(6) are given by CiF(C 2 p _1 ) ~ CiC|F(p _1 ) with constants C\ and C 2 . Knowing that outage 
event occurs when the instantaneous SNR pz falls below certain threshold, and since F(C 2 p~ 1 ) is 
the probability that pz falls below C 2 , the asymptotic error rate expression CiF(C 2 p~ 1 ) represents 
a scaled version of the outage probability. Furthermore, C\ and C 2 depend only on the system 
specifications (a, (5, 9 t , g\(9) and g 2 (9)) and the variation exponent d of F(z). As a result, we 
only need to know d in addition to F(z) to obtain the asymptotic error rate and do not need 
to express F(z) or the corresponding PDF f(z) in series expansion form. For simple practical 
channels like Nakagami-m, d = m can be seen by inspection of F(z) or f(z). For many channel 
distributions, the diversity order can also be obtained by solving lmx^o F(rz)/F(z) = r d or 
lim^o f(rz)/f(z) = t^ 1 using L'Hopital's rule, or approximated numerically by evaluating 
lim^ log F(z)/ log z = d or lirn^o log f(z)/ \ogz = d — 1. Therefore, the asymptotic error 
rate is related to the channel distribution in many useful ways which complement the PDF-based 
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approach in [2]. In Figures 1 and 2, we compare the BERs of DPSK and BPSK under Nakagami- 
m fading obtained through Monte Carlo simulation, with their approximations given by (3) and 
(6). We observe that the results given by (3) and (6) match their corresponding simulation results 
within 0.1 dB at error rates of 1CT 6 . In addition, it can be seen that our approach gives better 
approximation to the Monte Carlo simulation results than [2], most noticeably for moderate 
values of average SNR. 

There are some practical cases of instantaneous error rate P e (pz), which do not fit ASlc 
but can be expressed as linear combinations of exponential-mixture functions of pz in the 
form J * g 2 (9) exp(— pz/gi(9))d9 given in ASlc, such as the instantaneous BER of Gray- 
coded M-PSK [15, Section 8.1.1.3]. There are also practical P e (pz) which fits ASlc, but 
can be expressed as a linear combination of exponential-mixture functions, each with a much 
simpler #2(0) function than that of P c (pz) itself, such as the instantaneous SER of square M- 
QAM. For these cases, asymptotic characterization of P c (p) through linear combination becomes 
necessary. Specifically, for an instantaneous error rate given by P e {pz) = Ylj=i a j^e,j(pz) 
where {%}/ =1 are constants and {P cj (pz)}j =1 are (simpler) expressions satisfying ASlb or 
ASlc, we first determine E[P c j(pz)] ~ CijF(C 2 jP~ 1 ) using the methods in this section. If 
{%}/=! are all positive, we have P c (p) ~ J2j=i a jCijF(C2jP~ 1 ); if any of {%}/ =1 is negative, 
Pc(p) ~ J2j=i a jCijC2jF(p~ 1 ) can be established. We omit the derivations due to lack of space. 

VI. Asymptotic Error Rate Performance under Diversity Combining 

In this section, we establish an extension and application of the results derived in Section V-B, 
by analyzing the asymptotic error rate performance at high average SNR for several diversity 
combining schemes. The results are especially useful when the diversity branches have non- 
identical fading distributions and when the average error rate expression is not available in 
closed form. 

We consider a system in which the receiver has the channel state information (CSI), and 
employs N diversity branches with independent but not necessarily identical fading distributions. 
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Particularly, the n-th branch is assumed to have a channel power gain z n with CDF F n (z) = 
z dn l n (z) where l n (z) is slowly varying at 0. We derive asymptotic expressions for E[/3e~ apZc ] and 
EI/q 6 ** 92(6) exp(— pz/gi(6))d6] for large p in terms of the system specifications (a, (3, 6 t , g\{9) 
and <?2(#)X P> {d n }n=i an d On(')}^=i' wnere z c is the channel power gain after combining. In 
order to address this problem with a unified approach, we first determine the asymptotic CDF 
of z c near the origin in the form F c (z) ~ z dc l c (z) with l c (z) being slowly varying at for each 
specific diversity combining scheme by expressing d c and l c (-) in terms of {d n }^ =1 and {l n (-)}n=i 
respectively. Moreover, we express the CDF of the combined channel F c (z) in terms of {d n }^ =1 
and {F n (-)}£Li- This will lead to characterizing E[f3e~ apZc ] and E[j^ g 2 {6) ex-p(-pz/ gi (6))d6} 
in terms of the system specifications, p, d c and l c (-) (and alternatively F c (-)) using the same 
method as in Section V-B. 

A. Maximum Ratio Combining (MRC) 

For MRC with independent branches as well as a number of cooperative relay systems [19], 
the channel power gain can be expressed as the sum of independent random variables: z c = 
J2n=i z n- Define C n (s) := / °° e~ sz dF n (z) and C c (s) := J °° e~ sz dF c (z) to be the Laplace- 
Stieltjes transforms of the respective distributions. Based on the convolution property of Laplace 
transform, we have C c (s) = Yln=i £n{s). On the other hand, it follows from Tauberian theorem 
(Theorem 2) that C n (s) ~ T(d n + l)s~ d "/„(s _1 ) as s ->■ 00, and therefore C c (s) ~ f[n=i r (^™ + 
1) . s-EiLi*. YlLiU 8 ' 1 )- K is easy to verify that Un=i T ( d n + 1) • Un=i U^ 1 ) is slowly 
varying at 00 as a function of s. It follows from Tauberian theorem (Theorem 2) that 



near the origin, where we have simplified z^=i d ^ FJ^ =1 l n (z) as Y\n=i ^n(z), and the first factor 
on the right hand side is due to the variation exponent — J2n=i of C c (s) at 00. We will relate 
the outage of the combined channel in (7) to the average error rates in Section VI-D. 
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B. Equal Gain Combining (EGC) 

For EGC with independent branches we have z c = (J2n=i \^n) 2 /^- Define y n = ^fz^ and 
y c = J2n=i Vn- It is easy to derive that y n has a CDF given by G n (y) = y 2dn l{y 2 ), and show that 
l(y 2 ) is slowly varying as a function of y at 0. Consequently, by using the same method as in 
Section VI-A, we can derive the asymptotic CDF G c (y) ~ [r(2 En=i d n + I)] -1 Un=i T ( 2d n + 
1) • y 2 ^ =1 dn Y\n=i ln(y 2 ) of y c near y c = 0, then using z c = y 2 /N, we get the asymptotic CDF 



F c (z) 



N 



r 2^4 + 1 



n=l 



-1 



N N 



l[T(2d n + l)-l[F n (Nz) (8) 

n=l n=l 

of z c near the origin, which will be related to the asymptotic error rates in Section VI-D. 

C. Selection Combining (SC) 

For SC with independent branches we have z c = max ri=12i ... i Ar{-2„}, then it follows that 

N N 

F c (z) = l[F n (z) = z ^n=i^Y[l n (z), (9) 

n=l n=l 

which will be used to derive the asymptotic error rates next. 

D. Asymptotic Error Rate Expressions 

For all three cases analyzed above, we can verify that F c (z) is regularly varying at with 
exponent d c = J2n=i d 'i > 0, since the same holds for F n (z) with exponent d n > 0. In Section 
V-B we established the asymptotic average error rates in terms of the CDF F(z) in closed- 
form. We can directly apply these results here to the combined CDFs of the respective diversity 
combining schemes, and obtain 

E[f3e~ a ^] ~ f3T(d c + l)F c {a' l p- 1 ) (10) 
J\ 2 (6)exp(--^)d9 ~ (r(d c + l) £ g 2 (9)gl°(9)d9^ Fcip' 1 ) (11) 



E 



as p — > oo, where d c is the same for all three combining schemes, and F c (-) is given by (7), 
(8) or (9) for MRC, EGC and SC respectively. We have thus established a unified approach 
to evaluate the asymptotic error rate at high average SNR for MRC, EGC and SC, given the 
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conditions that the channel power gain of each branch has a CDF which is regularly varying 
at and the instantaneous error rate can be expressed as linear combination of exponential or 
exponential-mixture functions. Figures 3 and 4 show the BERs of BPSK under different diversity 
combining schemes obtained through Monte Carlo simulation, as well as their approximations 
obtained through (11) together with (7), (8) and (9). The accuracy of our approach is corroborated 
by the closeness (within 0.2 dB at error rate of 10~ 6 ) between the simulation results and 
analytical asymptotic approximations. Like observed in Figures 1 and 2, our approach gives 
better approximation than [2] to the simulation results, especially when the average SNR is not 
significantly high. 

In addition to the simple relation between the asymptotic error rate and the channel dis- 
tribution, it can be seen that our approach based on Tauberian theorem enables analysis of 
the performance of a communication system involving multiple channels, in which the regular 
variation property of the distribution of the overall effective channel is inherited from the channel 
distributions corresponding to the constituent parts of the system. Also, the asymptotic CDFs 
given by (7) and (8) can make high-SNR substitutes of the numerical inversion method in [20] 
to compute outages in fading channels, as long as different diversity branches are independent. 

VII. Conclusions 

In this paper, we investigate the relationship between the high-SNR asymptotic average error 
rate and outage by establishing that the instantaneous error rate being upper bounded by an 
exponential function of the instantaneous SNR is sufficient for outage events to dominate error 
rate performance. We do this by establishing the regular/rapid variation of the CDF of the 
channel power gain with exponent d at to be a necessary and sufficient condition for the error 
rate to exhibit a diversity order of d. This implies the equivalence between the error-rate-based 
diversity order and the outage-based diversity order, thereby characterizing the conditions under 
which the outage event dominates high-SNR error rate performance. For the case of finite d, 
the constant SNR offset between different communication systems over the same fading channel 
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with exponentially bounded instantaneous error rates is revealed. For instantaneous error rates 
given by a mixture of exponential functions of the instantaneous SNR (ASlc), we express the 
high-SNR asymptotic error rates directly as multiples of the outage probability (channel CDF for 
small arguments). Furthermore, we derive the asymptotic error rate for instantaneous error rate 
satisfying ASlc, considering different diversity combining schemes with the CDF of the channel 
power gain of each diversity branch being regularly varying at 0. All these results exhibit the 
convenience and effectiveness of Tauberian theorem as a tool to analyze the asymptotic error rate 
performance related to diversity under fading, thus conveniently generalizing and complementing 
the PDF-based approach in [2]. Numerical results show that our approach gives more accurate 
approximations than [2]. 



We prove this for H(x) being slowly varying at oo only, since the case of slow variation 
at can be obtained by considering H(x^ 1 ). By definition we have lim^oo H(tx)/H(x) = 1 
for t > 0, which implies that for any e 1 > 0, there exists sufficiently large si such that for 
all x > si, \H(tx)/H(x) — 1| < €\. Without loss of generality we assume t > 1, and since 
< hi < H(x) < h 2 < oo for sufficiently large x (say x > s 2 ), it follows that for any 
e 2 = eih 2 , there exists sufficiently large s 3 = max{si, s 2 } such that for all x\ > s 3 and x 2 > s 3 , 
\H(x 2 ) — H(xi)\ < e 2 . This follows from multiplying \H(tx)/H(x) — 1| < ei by H(x). Therefore 
the Cauchy criteria for the existence of limit is satisfied, and H(x) should converge to a constant 
in [hi, h 2 ] as x — > oo. 



We use Theorem 1 to prove the relation between (2b) and (2c), starting from the claim that 
(2b) implies (2c). For < d < oo, (2b) implies that P c (p) has a variation exponent — d at oo, 
and thus can be represented as 



for some a > 0, c(p) — > c e (0, oo) and e(p) — > — d as p — > oo, using Theorem 1. Therefore 



Appendix A. Proof of Proposition 2 



Appendix B. Proof of Proposition 3 




(12) 
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lim iogPe(P) = lim \ogc( P )+ir-fdu = iim jr-^du 

p^oo \ogp P^oo logp P^°o logp ^) 

= l im !kl( p -i)-i = i im e ( p ) = 

p— >oo p p—yoo 

where in the second equality we used e(p) — y —d, and in the third equality L'Hopital's rule. We 
have thus shown that (2b) implies (2c). 

We next prove that (2c) implies (2b) given the mild additional assumption that the limit in 
(2b) exists. The existence of the limit in (2b) implies that lim^oo P c (rp)/P c (p) = r mi for some 
mi E [— oo, oo] [8, Lemma 1, p. 275]. We would like to show that mi = —d. Clearly, P e (p) can 
be represented by (12) for some a > 0, c(p) — y c G (0, oo) and e(p) — y mi as p — y oo. Similar 
to (13), it can be derived that lim^oo (log P e (p))/ (logp) = mi. Given the condition that (2c) 
holds, it follows that mi = —d, and thus (2b) also holds. 

To show that (2a) is a special case of (2b), from (2a) we get P e (Tp) = G~ d T~ d p~ d + o(p~ d ), 
then 

lim em = Um G ;y + fp-'y'-' = T -' 

p^°° P c (p) p^°° G- d + o(p- d ) I p- d 
i.e. (2b) holds. An example for why (2b) does not imply (2a) can be seen for the case P c (p) ~ 
ci(log p) C2 p~ d with ci,c 2 > 0, as brought up in [6], [7]. Hence (2b) is more general than (2a). 

Appendix C. Proof of Theorem 3 

We begin with assuming the regular/rapid variation of F(z) with exponent d at 0, and proving 
the regular/rapid variation of P e (p) with exponent — d at oo, by analyzing the upper and lower 
bounds of P e (p). Due to ASla we have 

/•oo roo 

P c (p)=/ P c (pz)dF(z)< l3e- apz dF{z) = Y c Mp) (15) 
Jo Jo 

and since P e (-) is non-negative and monotonically decreasing, 

p/p rv/p 
P e (p) > / P c (pz)dF(z) > / P^dFiz) = P C ( V )F(V/P) = Pc,lb(p) (16) 
Jo Jo 

for any constant rj > 0. It follows from (15) and (16) that 
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log -r e,LB (p) < logPc(p) < logP c ,U B (p) 

logp ~ logp _ logp 

For d — oo, it can be easily seen that P c ,lb(p) has variation exponent — oo at p = oo and 
thus satisfies (2c). Also, based on Tauberian theorem (Theorem 2) with the correspondences 
F(z) = H(x) and ap = s, P e ,uB(p) has variation exponent — oo at p = oo and thus satisfies 
(2c). By taking the limit of (17) as p — > oo and using Squeeze theorem (a well-known theorem 
stating that if \im x ^ a H 1 (x) = lim x ^, H 3 (x) = h and H 1 (x) < H 2 (x) < H 3 (x) in some 
neighborhood of a, then \im x ^ a H 2 (x) = h), it follows that P c (p) should satisfy (2c). 

For < d < oo, let F(z) = z d l(z) with l(z) being slowly varying at 0. Based on (15) and 
Tauberian theorem (Theorem 2) we have 

P e ,u B (p) ~ PT(d+ l)(ap)- d Z(a-V _1 ). (18) 
Also, it follows from (16) that 

p c Mp) = v d p- d P c (v)Kv/p)- (19) 

It is easy to verify that both P Cj ub(p) and P e ,LB(p) satisfy (2b), thus also satisfy (2c) due to 
Proposition 3, then P e (p) should satisfy (2c) due to (17). This completes the sufficiency part. 

We now show P e (p) being regularly/rapidly varying with exponent — d at oo implies F(z) 
being regularly/rapidly varying with exponent d at provided that F(z) is a Karamata function 
(at 0). Based on [8, Lemma 1, p. 275], it can be confirmed that lim 2 ^ F(tz)/F(z) = r" 12 
with m 2 £ (0, oo] being the variation exponent of F(z). Like in the sufficiency part of the 
proof, it can be derived that both P e ,uB(p) and P e ,LB(p) have variation exponent — m 2 , then (17) 
becomes — m 2 < — d < — m 2 , leading to m 2 = d. Therefore F(z) is regularly/rapidly varying 
with exponent d at 0. 

Appendix D. Proof of Theorem 4 

We first derive the asymptotic expression for P e (p) assuming AS2b, F(z) = z d l(z). Theorem 
3 implies that P c (p) = p~ d r(p) where r(p) is slowly varying at oo. Let h(p) := P e (p) / 1 (p~~ d l(p~ v )) 
It is easy to show that h{p) is slowly varying at oo. On the other hand, from (15) and (16) 

DRAFT 



20 



we have P c ,lb(p)/(p^(p" 1 )) < h(p) < P c ,ub(p)/(p'^(p" 1 )); also, based on (18) and (19) as 
well as the definition of slow variation, we have linip^oo P ei uB(p)/(p^ d ^(p -1 )) = /3a~ d T(d + 1) 
and lim^oo P etLB (p) / (p^ d l(p^ 1 )) = P c (^)^ d respectively. It follows that for sufficiently large p, 
h(p) is bounded between two finite positive constants. Based on Proposition 2, h(p) converges 
to a finite positive constant as p — > oo, say h(p) ~ c\ G (0, oo), and hence P e (p) ~ c iP d Kp 1 )- 
Consider now the derivation for the asymptotic expression of F(z) assuming P c (p) = p~ d r(p) 
and that lim 2 _^ F{tz)/F{z) exists for r > 0. It follows directly from Theorem 3 that F(z) 
must be in the form F(z) = z d l(z) with l(z) slowly varying at 0. Based on the sufficiency 
part of the proof (i.e. the asymptotic expression of P c (p) given F(z) = z d l(z)), we must have 
P e (p)/(p~ d l(p~ 1 )) = r(p)/l(p^ 1 ) ~ ci as p — > oo, with constant ci G (0,oo). Consequently, 
l(z)/r(z~ v ) ~ cf 1 as z — > 0, and F(z) ~ c 2 z d r(z~ 1 ) where c 2 = cf 1 G (0, oo) is a constant. 
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Fig. 1. BERs of BPSK and DPSK under Nakagami-m fading with m = 2 
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-* — BPSK by simulation 
i DPSK by simulation 
-A — BPSK by approximation 

BPSK by results in [2] 
-e — DPSK by approximation 
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Fig. 2. BERs of BPSK and DPSK under Nakagami-m fading with m = 3 
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Fig. 3. BERs of BPSK under MRC and SC involving three Nakagami-m fading branches with m = 0.5, 1 and 1.5 
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Fig. 4. BERs of BPSK under EGC involving three Nakagami-m fading branches with m = 0.5, 1 and 1.5 
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